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Abstract 

We investigate some structure scalars developed through Riemann 
tensor for self-gravitating cylindrically symmetric charged dissipative 
anisotropic fluid. We show that these scalars are directly related to 
the fundamental properties of the fluid. We formulate dynamical- 
transport equation as well as the mass function by including charge 
which are then expressed in terms of structure scalars. The effects 
of electric charge are investigated in the structure and evolution of 
compact objects. Finally, we show that all possible solutions of the 
field equations can be written in terms of these scalars. 
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1 Introduction 

It is believed that anisotropy plays a vital role for understanding the gravi- 
tation of those objects which have higher densities than neutron stars. Many 
phenomena such as the solid core, phase transition, mixture of two fluids, slow 
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rotation and pion condensation can generate anisotropy in the star models 
PP. Lemaitre [2] was the first who gave the idea about the tangential and 
radial pressures. Since the pioneering work of Bowers and Liang [3], there 
has been extensive literature devoted to the study of general anisotropic rel- 
ativistic configuration both analytically and numerically [I]. 

Recent literature indicates interesting consequences of the inclusion of 
an electromagnetic field to discuss gravitational collapse. Bekenstein [5] 
extended the Oppenheimer-Volkoff equations of hydrostatic equilibrium [6] 
from the neutral to the charged case. Herrera et al. [7] provided a set of equa- 
tions for the physical interpretation of models for collapsing charged spheres. 
Nath et al. [8] discussed charged gravitational collapse and concluded that 
electromagnetic field increases the formation of naked singularity. Sharif and 
his collaborators [9] investigated the effects of electromagnetic field on dif- 
ferent aspects of spherically/cylindically and plane symmetric gravitational 
collapse. 

The orthogonal splitting of the Riemann tensor was first considered by 
Bell [ID] . Herrera and his collaborators [IT] followed this idea to develop a 
relationship between structure scalars and the fluid properties. Also, they 
[T2] analyzed the structure scalars for charged dissipative spherical fluids. 
Structure scalars, Xt, Xtf, Yt, Ytf have important properties such as Xt 
is the energy density of the fluid, Xtf controls inhomogeneity in the fluid, 
Ytf describes the effect of the local anisotropy of pressure as well as density 
inhomogeneity of the Tolman mass and Yt turns out to be proportional to 
the Tolman mass density for systems in equilibrium or quasi-equilibrium. 
This ultimately provides a relationship between structure scalars and energy 
density inhomogeneity. 

In a recent paper, Herrera et al. [13] discussed cylindrically symmetric 
relativistic dissipative fluids based on structure scalars. Here we take the 
effect of electromagnetic field to study the structure scalars with the same 
configuration. The paper is organized as follows. In the next section, we 
describe the Einstein-Maxwell field equations, kinematics and the Weyl ten- 
sor. Section 3 is devoted for the structure scalars and conservation laws. In 
section 4, we formulate the dynamical equations in terms of mass function 
and couple with transport equation. We discuss the possible static charged 
cylindrically symmetric solutions for the anisotropic fluid in section 5. In the 
last section, we summarize the results. 
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2 Charged Anisotropic Dissipative Fluid Cylin- 
ders 

Here we review the basic general equations and some definitions. The general 
cylindrically symmetric spacetime is given as 

ds 2 = ~A 2 (t, r)(dt 2 - dr 2 ) + B 2 (t, r)dz 2 + C 2 (t, r)d(P 2 , (1) 

where — oo ^t^oo, O^r, — oo<z<oo, 0^0^ 2ty. We assume 
that the collapsing cylinder is filled with anisotropic and dissipative fluid for 
which the energy-momentum tensor is 

T$ = (n + P r )V a V p + P r g aP + q a Vp + qpV a + U a p, (2) 

where U a p = (P z - P r )S a Sp + (P^ - P r )K a K/3, q a = qL a , // is the energy 
density, P r , P z , P^ are the pressure in the radial, z and <j) directions, respec- 
tively, q a is the radial heat flux and V a is the four velocity Also, S a , K a 
and L a are the unit four-vectors with 

V a V a = -1, L a L a — S a S a — K a K a — 1, 
V a L a = V a S a = V a K a = S a K a = 0. 

For comoving coordinate system, we have 

V a = -AC L a = ASl S a = B6 2 a , K a = Cd 3 a . 

The energy-momentum tensor for electromagnetic field is 

t7 ] = i - l^TsgJ) , (3) 

where F a/3 = <pp j0l — <f> a> p is the Maxwell field tensor and <p a represents the 
four potential. The Maxwell field equations are 

Ff = ^J a , F [aft7] = 0, (4) 

where — 47r is the magnetic permeability and J a is the four current. In 
comoving coordinate system, we can assume that the charge per unit length 
of the system is at rest, so the magnetic field will be zero. Thus one can 
choose the four potential and four current as 

<P a = <P5l J a = (V a , 
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here (p is the scalar potential and ( is the charge density, both are functions 
of t and r. The only non-zero component of the Maxwell field tensor is 



F ( 



01 — -^01 



where prime is the differentiation with respect to r. Using these values, the 
Maxwell field equations become 



+ *[b + c—a) = *^ a - (5) 

d f 1 d<f>\ ( I d(f)\ ( B C 2A\ n 

m{T*fr) + l v fi + c + Tj =0 ' (6) 

where dot means differentiation with respect to t. Integration of Eq. (j5]) yields 

2s A 2 



BC ' 

where 

FT 

s( r ) = 2tt / C,ABCdr, (7) 
Jo 

is the total amount of charge per unit length of the cylinder found through 
the conservation equation, = 0. Also, </>' identically satisfies Eq.flHJ). The 
Einstein-Maxwell field equations yield 

, 2 / s 2 \ ( B C\ A BC B" C" 

kA 2 fi + 



2ttB 2 C 2 J \B C A BC B C 

A' ( B' C'\ B'C 
+ T 75" + 77 - 7777' ( 8 ) 



-KqA 2 



+ T U + 7? . ( 9 ) 



, s 2 \ _B_C AfB C\_BC 

1 r 2tvB 2 C 2 ) ~ B C + a[b + C BC 



A' 


IS 


C'\ 


B' C 


a' 




-c) 


+ B C 



+ T 77 + 77 ( 10 ) 
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There are four kinematical variables for the description of fluid, i.e., ex- 
pansion, acceleration, shear and rotation. Since we are using irrotational 
fluid, so the rotation will be zero. Rest are defined as 

6 = V£, a a = V a .pV fi , a a/3 = V( a 0) + a {a Vfs) - ~Bh a/3 , 

where h a p = g a p + V a V/3. Using Eq.fll]), these quantities turn out to be 

^ 1 [A B d\ A' 2 Q (A'V 

6 = -l- + - + -l, = — , a=aa «=\^), a a = aL a . 

The shear tensor can also be expressed as 

C a j3 = °"s (s a S/3 - + a k (^KaKp - -h aj3 

= ^(a 2 s -a s a k + a 2 k ), (13) 



where 



1 (B A\ l(C A 

as -A\B-A > ak -A\C~A 



The Weyl tensor is defined as 



which may be decomposed in its electric and magnetic parts as 

E aP = C^sV'V 5 , H a p = ^n avep Cf 5 V v V s , (14) 

respectively, where 7] auep is the Levi-Civita tensor. These can also be written 

as 

1 1 

E aP = E s (S a Sp--h a p) + E k (K a Kp--h al3 ) 1 H aj3 = H(S a K^ + SpK a ), (15) 

where 

tp - 1 n 1 r tp - 1 r> 1 r u — C ° 313 
~ A^EP 202 _ A 2 101 ' ~~ ^2(32 0303 ~ 342 Uolol » ^ ~ ^A 2 ^' 

; (i6) 

The components of the Weyl tensor C0202, C0101, Co303 3 C0313 are given in a 
recent paper [13]. 



3 Structure Scalars 

In this section, we formulate structure scalars for the charged fluid from 
the orthogonal splitting of the Riemann tensor [13] . For this purpose, the 
following tensors are defined 

where R^^s = ^VepysR^- These can be expressed in trace and trace free 
parts as 

Y a /3 = -Yrha/3 + Y s (S a Si3 — —h a p) + Y k (K a Kp — ~h a p), 

X a /3 = -Xxha/3 + X s (S a Sp — -h a p) + X k (K a Kp — -h a p), (17) 

Using the field equation (l8l-(fT2l) with (fl6l) . we obtain Y a p and X a p in terms 
of physical variables 

Yt = ^(V + Pz + Pt + Pj + J^, X T = Kfi+^^, (18) 

Y s = "(P,_p r )_-i£-, Y k = E h -%(Pt-P r ) |S 



2 y ' B 2 C 2 '' 2 B 2 C 2 ^ 

k As 2 k As 2 

X R = —En — —\Pz — Pr) „r, , Xk = —E k — —(Peh — Pr) 



2 - - rj B 2 C 2 '' 2 B 2 C 2 ' 

(19) 
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The conservation law, T°f = 0, yields 

(A* + Q(ji + P r ) + + a a q a + a a ^ + ^QU a a = 0, 

4 ss' 
h^{P r0 + + <£) + (jjl + P r )a a + -Qq a + ay - = 0. (20) 

The last equation can be written in an alternative form as 



P X + Q* - — 
r q A 



+ (n + P r )a 



1 ss' 
~ -^ + ^-^-^AB^ = ' (21) 

where /* = f, a V a , f' = f, a L a . There are two important differential equa- 
tions relating the Weyl tensor to different physical variables. Herrera [13] 
found these relations for cylindrically symmetric spacetime which are gener- 
alized to charge distribution as 

-(Y s + Y k -X s - X k ) x - 3(Y S -X s )^-- 3(y fe -X k ) C 



AB v " *' AC 
-6H(a s - a k ) = n(2fi + P r + P z + + 3«(^ + P r )a 

+ 2 Kg (0 -a s - a k ) + 3nq* + ^ 2(J2 - ^ - s^j , (22) 
(2Y S -Y k - 2X S + X k y + 3(Y S - X s )^ + 3a(Y s -Y k -X s 



+X k ) + 6H(Q - a k ) + 6H* = -k{h - P r - P z + 2P^ 

(.f-/). (23) 



C 3ks 
-3k(P - P r )-r^ + nq{G - a s + 2a fc ) + 



AC y 7rAB 2 C 2 

4 Mass Function and Dynamical- Transport 
Equation 

Here, we develop equations that govern the dynamics of non-adiabatic cylin- 
drically symmetric collapsing process. For this purpose, we define the velocity 
U = % = C*. Using Eq.(HnP, we get 



C C B AB BC B'C A'B'\ ks 2 
U* = a— - nP r C I — 



A r A 2 \ B AB BC BC AB 2<kB 2 C 
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which turns out to be 
IT 



0202 



R 



A 2 



2323 
~C 2 ~ 



+ 



KS 



2nB 2 C 



Solving it for a and substituting into Eq.(l2TT). we obtain 



(fX + P r )U* = -(fi + Pr) 



:p, ~ Pr 



KP r C - 



C_ ( R 

B 2 



H)202 

a 



R 



2323 



AB ^ Pr > AC 



+ a fe -46)g] + 



ss'C 
tiB 2 C 2 



(fi + P r 



C 2 

91 
~A 

2nB 2 C 



~A 



[Pi 



[cr. 



(24) 



This equation yields the effect of different forces on the collapsing process 
and have the "Newtonian" form as 

Force = Mass Density x Acceleration. 

The term on the left is the density (the inertial or passive gravitational mass 
density) multiplied by the proper time derivative of the velocity U . The terms 
on the right represent the force which is the contribution of four forces: the 
gravitational force (the pressure gradient plus the anisotropic contribution), 
the contribution from the dissipation and the contribution from the charge 
term. 

Now we define a mass function m similar to the spherically symmetric 
case. For this purpose, we write the term R ™ 2 
scalars as 



R ^i 3 in the form of structure 



R 



0202 



R 



2323 



-{Y T 
3 V 



A 2 C 2 
Using Eqs.flU]) and (02}, it follow that 

Rmno KB 2 



X T + X s + 2Y S + X k -Y k 



M3202 



12323 
C 2 



| + 2P r 



Pa, 



B 2 



P z + ^)+^r(E s -2E, 



As 2 

w < 25 > 



Inserting this value in the first square brackets on the right side of Eq. (I24j) . 
we obtain 



C (R 



B 2 



t-0202 

~A^~ 



R 



2323 



C 2 



kP t C 



K 



+ (^-P^-P r + 2P Z )C 



6 

- f <* 



2E k ) + 



As 2 

eJc' 



(26) 
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To define a mass function, we have followed the procedure of [7J[T3]- We 
have compared Eg. (126 p with the corresponding equation in the spherically 
symmetric case [TJ. We have also assumed that the pressure effect remains 
the same as in spherically symmetric case which provides the definition of 
mass function as a possible extension of the Misner-Sharp mass function to 
the cylindrically symmetric case given by 

m = ^(/i -P+- P r + 2P Z ) - ^(E s - 2E k ) + (27) 

3 B z 

Making use of Eq.( 1T9"j) . this can be written as 

= f 0" + P* - 2Pr + P.) ~ (Ys ~ 2Y k ) + ig, (28) 

or 

— = -(//- 3P^ + 3P Z ) - (2X k - X s ) + — . (29) 

The mass function in terms of electric charge can be found by using 
Eqs.flD and flUD-flU as follows 

(3Y S - 3Y k + X s + X k y = nrf + nq{a k - - 2a s ) - ^-{Y s + X s ) 

AB 

3C' 

■ J (Y k - X k ) - QH(Q - a,) - QH* - 3a(Y s - Y k - X s + X k ) 



AC 

3ks ( , sC 



+ ^aWc> {'--)■ < 30) 

Adding Eqs.(J2SJ and ([29]), it follows that 

( KfI + 3Yk -3Y s -X k -X; 



Applying the operator f (/t = f\ a L a ) on both sides of the above equation 
and then substituting in Eq.( )30|) . after some manipulation, we obtain 

/6mV B' 

i-^j = 3(Y s + X s )— + 6H(Q-a s )-3(Y k -X k ) 
C 

x — + Kq(2a s -a k + Q) + 6H* + 3a(Y s - Y k 

v , ks f , sC AsB' 
s+ k) nAB*C*V C B 
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Integration leads to 



rn 



6 



3(Y S + X S )^- + 6H(Q - a s ) - 3(Y k - X k ' ( 



AB 



+ nq(2a s - a k + 9) + 6H* + 3a(Y s -Y k -X s + X k 



AB 2 



TV 



sC AsB' 



C 



B 



dr + 



C 3 7(t) 



AC 

KS 

C 2 



6 



(31) 



where 7 is an arbitrary integration function of t. This mass function shows 
its dependence on different factors, in particular, on electric charge. Inserting 
the values of X s , Y s , X k , Y k from Eq. (1191) . it can be expressed in terms of 
physical variables and the Weyl tensor 



m 



C 3 
2 



.4 



n(p r - p.)^= - 2 -^- + 2H(e -* s ) + ^- 

KS 



AB AC 

C 3 f 

x (2a s -a k + Q) + 2a{E s - E k ) +2H*} dr - — 

C 3 7 (t) 



TV 



B 2 C 2 



sC AsB' 



s — 



C 



B 



dr + 



6 



(32) 



This gives the contribution of anisotropy, Weyl tensor, electric charge and 
dissipation. 

The transport equation for dissipative fluids is given by [1] 



rh^Vq^ + q a = -Kh ap (Tp + Tap) - -KT' 



(kt 2 J . 



(33) 



Here K, T and r indicate thermal conductivity, temperature and relaxation 
time, respectively. The only one independent component is 



rq* + q = -K{T* + To) - -KT 2 q (-^) * - -qrB. 
Substituting Eqs.(f26D, (TJTD and flMD in (TjM]), [t follows that 



(34) 



(jl + Pr) 
1 - 



1 - 



KT 



Tfa + P r 



U* = -(» + P r 



P r C 3 



+ m 



1 

C 2 



KT 



t{h + P r 



~A 



-Pl + (P z -P r )^ + {P <j> -P; ( ' 



AB 



AC 
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a 

+ A 



KT 2 q f t \* q KT^ 1 . 5 ^ 



2r 



+ ^ + 2^c^ + P ^ [ l - T&TPrj ) • (35) 

This shows that gravitational attraction and electric charge on any fluid 
element will decrease by the same factor as the inertial mass density. 
We see that the charge increases the gravitational mass only if 

✓ >«Cfc + *)(-l + ;5 £y, (36) 

otherwise, it will decrease. This increase of gravitational mass causes rapid 
collapse. Bekenstein [5] noticed this strange effect for the Oppenheimer- 
Volkoff equations of hydrostatic equilibrium. We see that the charge does 
not enter into the term 1 .^v 7 ! . but it affects the gravitational mass and 

T(fl+P r ) & 

shows how thermal effects reduce the effective inertial mass. We observe that 
as T (ii+p ) — * 1' ^ ne i ner ti a l mass density of the fluid element tends to zero 



To] . This shows that there is no inertial force and matter would experience 

KT 



a gravitational attraction which causes the collapse. For < T , K + P \ < 1 



the inertial mass density goes on decreasing while T £+ P \ > 1 indicates the 
increase of inertial mass density. By the equivalence principle, there should 
occur increase or decrease of mass. Thus one can easily distinguish between 
expanding and collapsing mechanism during the dynamics of dissipative pro- 
cess. 

Assume that the collapsing cylinder evolves in such a way that the value 
of T £+ P \ increases and approaches to 1 for some region. During this process, 
the gravitational force term decreases and leads to a change of the sign of 
the right hand side of Eq. (l35]) . This would happen for small values of the 
effective inertial mass density and implies a strong bouncing of that part 
of the cylinder |16j . This phenomenon causes the loss of energy from the 
system and hence the collapsing cylinder with non-adiabatic source leads to 

the emission of the gravitational radiations. Notice that the term 1 , ^ T D t 

(in the inertial mass, gravitational force and electric charge) is related to the 
left of Eq.(153]). 
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5 Static Charged Anisotropic Cylinders 

Here we discuss all possible solutions of the field equations of the static 
charged anisotropic cylinders. The corresponding field equations are 



KA2 ( /i+ 2^) 



kA 1 R 



s 2 



2nB 2 C 2 



B" 




C" 


A' 


(B 


c\ 


b' a 


~~B 




C + 


A 




c) 


B c 




B' 






b' a 






-A 


B 


+ c 


h 


B C" 







2irB 2 C 2 J A C \A 

s 2 \ A" B" (A 1 

KA [P 6 + - -—: = — + 



2ttB 2 C 2 J A B V A / 

These equations, in terms of auxiliary variables, u = £ = V — %~-. 
can be written as 



s 2 



kA ^ + 2^B 2 ~C 2 ~) = - T /- T f + uZ + uV-Zv, (37) 



KA \ Pr ~2^C 2 ~) = "Z + UV + ZV, (38) 
^( P - + 2tW<p) = ^' + ^ (39) 



-,2 



Adding all these equations, we have 

u' + u^ + ur] = Y T A 2 . (41) 
From Eqs.fl38D-flIiD, it follows that 

u' + i' + e-^-^-iV = k(P^-P t .)A 2 + ^^, (42) 

uj' + v' + ^-uZ-ujv-^ = K (P z -P r )A 2 + ^^, (43) 

i' + e-ri' + rf = k{P*-P,). (44) 
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In terms of auxiliary variables, the scalars E s and become 
E s = — [-J + rf + rf + u£ - ujt] - £r)] , 

E k = ^[-u' + £ + e-^ + u V -t v ]. (45) 

Equations §£5$, yield 

Y S A 2 = -J + u£, Y k A 2 = -J + lot]. (46) 
Integration of the first equation yields 

where a is an integration constant, which means B = B(A) or £ = for 
any Yg. When we integrate second of the above equation, we obtain 

A = 7e V / , 

where 7 is another integration constant implying C = C(A) or 77 = rj(u) for 
any Yf.. Consequently we can express any (u),£, rj) in terms of each other. 
This leads to the similar result as that of paper [13] with the effect of charge. 
Thus any static anisotropic solution is determined by a triplet of scalars 
(Yfc, Y s , Xk) or (Yfe, Y s , X s ) as in the charged free case. Similarly, we can 
discuss the case of isotropic cylinders. 



6 Conclusion 

This paper investigates the effects of electromagnetic field on structure scalars 
of the cylindrically symmetric anisotropic dissipative fluid. The electric 
charge increases the inhomogeneity produced by local anisotropy. Following 
Herrera [13] . we have formulated dynamical as well as transport equations 
and also mass function for the charged cylindrical system. It turns out that 
the coupled dynamical-transport equation has an extra factor due to charge. 
If we take 7-^-5 ; = 1, the inertial mass and gravitational force vanish while 

T(fl+P r ) ' to 

the gravitational mass reduces for s' < skC(h + P r ) ^—1 + T £+ P ■ Further, 

we have discussed the static case in electromagnetic field which shows that 
any solution of the field equations can be expressed in terms of the scalar 
functions like charge free case. It is worth mentioning here that all our results 
reduce to charge free case p2] for s = 0. 
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